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Abstract 

There are two apparent puzzles connected with the two-body and three-body doubly charmed 
baryonic B decays. First, earlier calculations based on QCD sum rules or the diquark model 
predict B{B^ —> ~ —> BcN), while experimentally the former has a rate two orders 

of magnitude larger than the latter. Second, a naive estimate of the branching ratio 0(10“®) for 
the color-suppressed three-body decay B —> A^A~K, which is highly suppressed by phase space, 
is too small by five to six orders of magnitude compared to experiment. We show that the great 
suppression for the A^A~K production can be alleviated provided that there exists a narrow hidden 
charm bound state with a mass near the AcAc threshold. This new state that couples strongly to 
the charmed baryon pair can be searched for in B decays and in pp collisions by studying the 
mass spectrum of or AcAc- The doubly charmful decay B —> HcAc has a configuration 

more favorable than the singly charmful one such as B^ —> Acp since no hard gluon is needed to 
produce the energetic HcAc pair in the former decay, while two hard gluons are needed for the latter 
process. Assuming that a soft qq quark pair is produced through the a and vr meson exchanges in 
the configuration for B —> HcAc, it is found that its branching ratio is of order 10“^, in agreement 
with experiment. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 


Recently Belle has observed for the first time two-body and three-body doi^ly_charmed baryonic 

The measured 


Q,Ii- 


B decays in which two charmed baryons are produced in the final state 
branching ratios are 

B{B- ^ AtA~K-) = (6.5lo;^±0.8±3.4) X 10“^, 

B{B^ A+A-K^) = (7.9l|^ ± 1.2 ± 4.2) x 10“^, 
for three-body decays and 

B{B- ^ ^ H-7r+) = (4.81^;° ± 1.1 ± 1.2) x 10-^ 

H+A-)^(H+ ^ H-7r+7r+) = (9.3l|^ ± 1.9 ± 2.4) x IQ-^ 


( 1 ) 


B(B° 


( 2 ) 


for two-body decays. Taking the theoretical estimates (see e.g. Table III of |^), B{ 
B{Et- 




■ 1.3% and B{E'^ 
l^Tr^)/B(E'^ —> H“7r+7r^ 


-'tt^) ~ 3.9% together with the experimental measurement 


= 0.55 ±0.16 


it follows that 


B{B- 


°Ar) 


4.8 X 10 


-3 


B(^ 


■A“) Ri 1.2 X 10 


-3 


( 3 ) 


Therefore, the two-body doubly (±^med baryonic B decay B 
order 10“^, to be compared with 

=70 


BcB'^ has a branching ratio of 


B{B 


A+p) = (2.19l|];^^ ± 0.32 ± 0.57) X 10“® 


1.49 

B{B- Sc(2455)°p) = {3.67toit ± 0.36 ± 0.95) x 10“® 
for singly charmed baryonic B decays and [3,0| 

B{B° pp) < 2.7 X 10“^ 

BiB- Ap) < 4.9 X 10"^ 

for charmless baryonic B decays. Therefore, we have the pattern 


B{B^ AA) < 6.9 X 10 


-7 


( 4 ) 


( 5 ) 


BcB'^ (~ 10 -'^) > BcB (~ 10 “^) > B 1 B 2 (~ 10“0 


( 6 ) 


for two-body baryonic B decays. 

Using B(B^ —> A^p) as a benchmark, one will expect a branching ratio of order 10“^ for the 
charmless decay B —> B 1 B 2 after replacing the quark mixing angle Vcb by Vub, provided that the 
dynamical suppression for the latter is neglected. However, since the doubly charmed baryonic 
decay mode HcAc proceeds via h —> esc, while AeP via a 5 ^ cdu quark transition, the CKM mixing 
angles for them are the same in magnitude but opposite in sign. One may wonder why the BcB'^ 
mode has a rate two orders of magnitude l arg er than BcB. Indeed, earlier calculations based on 
QCD sum rules or the diquark model all predict that B{B —> HcAc) ~ B{B —> BcN), 
which is in violent disagreement with experiment. This implies that some important dynamical 
suppression effect for the BcN production with respect to EcAc is missing in previous studies. 

As for the three-body decay B —> AcAcK, its branching ratio is estimated to be of order 
10“®, which is extremely small due to the tiny phase space available for this decay and the color- 
suppression effect. The puzzle is that why the measured rate is much larger than the naive expec¬ 
tation? 
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FIG. 1: B~ —> A+A“iF“ as proceeding through (a) the internal VF-emission diagram, and (b) 
the dominant charmonium-like resonance Xcc- The blob in (b) shows where the strong decays take 
place. 

A crucial ingredient for understanding the baryonic B decays is the threshold or low mass en¬ 
hancement behavior of the baryon-pair invariant mass in the spectrum for B —> B 1 B 2 M: It sharply 
peaks at very low values. That is, the B meson is preferred to decay into a baryon-antibaryon pair 
with low invariant mass accompanied by a fast recoiled meson. Therefore, some three-body final 
states have rates larger than their two-body counterparts, e.g. ppK^ » pp, Apn^ ^ Ap, ScP'/r=*= S> 
This phenomenon can be understood in terms of the threshold effect, namely, the invariant 
mass of the dibaryon is preferred to be close to the threshold. The configuration of the two-body 
decay B B 1 B 2 is not favorable since its invariant mass is m^. In i? —> B 1 B 2 M decays, the 
effective mass of the baryon pair is reduced as the emitted meson can carry away a large amount of 
energies. The two-body decay pattern m also follows from the low-mass enhancement effect: The 
energy release is least for the B decay into two charmed baryons and becomes very large when the 
final-state baryons are charmless. 

Although the gross feature of the baryonic B decays can be qualitatively comprehended in terms 
of the near threshold effect, how to quantitatively evaluate their absolute decay rates and how to 
realize the low mass enhancement effect require detailed dynamical studies. In the present work we 
will focus on the doubly charmful baryonic B decays, namely, B —> HcAc and B —> AcAcK in Secs. 
II and III, respectively, aiming to resolve the aforementioned two puzzles connected with them. 
Sec. IV gives the conclusion. The evaluation of the delta functions occurring in the phase space 
integral is discussed in the Appendix. 


^ The three-body decay is usually referred to the nonresonant one. The relation Acpir ^ AcP is trivial as 
the former arises mostly from resonant contributions Q. 
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II. THREE-BODY DECAYS 


We consider the decay B —> K , which proceeds through the internal lY-emission dia¬ 

gram Fig. IHa). It turns out this diagram is factorizable. In the weak Hamiltonian approach, the 
factorizable amplitude reads 


A(R- ^ A+A-A-) = ^y,,K>2(A+A-|(cc)|0)(A-|(s-6)|R-), 


( 7 ) 


where {qiq 2 ) = ~ 75 )q 2 , and the effective Wilson coefficient 02 indicates that this decay is 

color suppressed. The matrix elements can be parametrized as 


(Ac(pi)Ac(p2)|(cc)| 0) = UA,(Pl) 


^ / 2x , •/2(g^) ( f 2^ , 53 ( 9 ^' 

mQ )7m + *Trz—- \9i[q + 


2mA, 


2mA, 


q^l 75 


{K-{pK)\{-sb)\B-{pB)) = F^^{q^){pB+PK)^+{Fi^{q^)-F^^{q^)) 


niu — m 


K 




qfi, 




( 8 ) 


with q = Pb — Pk = Pi + P2- la terms of the form factors, the decay amplitude has the expression 


with 


A{B ^ A+A^ A ) = ^VcsV*h a2UA, [a^K + b- {c^k + d)75]aA, 


a = 2Ff^iq^) Mq^) + f 2 {q^) 


b = 2F^^{q^)f2{q^){p2-pi) ■PK/{‘2mx, 
:^BK 


(9) 




d = 2mA,5i(9^) 


7 BK / 


(g^) + (Fo^^(g^)-Fr(g^)) 


BK/ 


— m 


K 


( 10 ) 


lll| where the 


+9z{q)Fo (g )(mB-m^)/(2mAj. 

There are numerous estimates of the B ^ K transition form factors. We will follow 
form factors are evaluated using the relativistic covariant light-front quark model. 

Due to the heavy mass of Ac, the phase space of the A+AjA“ decay is about a hundred times 
smaller than, say, that of the ApTr"*" jl2l| . The A^A“ form factors, if any, can therefore be taken as 
constants whose values are determined at the threshold over the phase space. To achieve a rate that 
is at least comparable with that of the ApTr"*" whose branching ratio is of order 3 x 10“® [l^ . one 
would need the A+ A.~ form factors to be more than one hundred times larger than those of the ApTr'*' 
near the Ap threshold, which is quite unlikely since the Ap form factors have their maximum values 
already about 0(1). Besides, O(IO^) form factors would just give a rate of A+A“A“ comparable 
to ~ 0(10“®), not to mention the remaining factor of O(IO^) difference between the rates of Apvr"’' 
and of A+A“A“. Therefore we conclude that the suppression from the A+A^A” phase space is so 
strong that A+A“ pair is unlikely to be produced dominantly through the direct three-body decay 
processes. The great suppression, however, seems to hint strongly that a cc-content resonance with 
the width comparable to the nearby resonances like ^(4415) could be located around the threshold 
of A+A- (~ 4.6 GeV), and the whole process takes place dominantly via the charmonium-like 
resonance as shown in Fig. Hb). 


4 













Let us assume the resonance, Xcc-, exists with a mass ^ 2mAc a width Tx^^- Let 

us further assume that this resonance is a spin-1 particle with = 1“ or I"*", as inspired from 
the observation that all the charmonia near the A^A“ threshold are spin-1 particles. The decay 
amplitude then reads 


AxAB- ^AtKK-) = ^VMa2{K-\s^^{l-^5)b\B-) 


^/2 


xmx,Jx,, 


( 


-g^^ + 






iuA,{pi)M^vj^^{p2), ( 11 ) 


where 


when Xcc is a vector {Xcc 


My 


hXXV 


Iv + 


2mA^ 


o-i/p , 


V), and 


My 


^Ac Acj 4. _ 


lu- 


2mA. 


75 , 


( 12 ) 


(13) 


when Xcc is an axial-vector {Xcc = A) particle, fx^c is the decay constant for Xcc, and 
and represent the dimensionless AcAcXcc strong couplings. Since the allowed phase space 

is very small, the strong couplings can effectively be treated as constants within this region. The 
decay constant fx^c comes from the factorization of the amplitude of B~ XccK~, followed by 
the strong decay Xcc —> A+A“. Since the chirality-flipping baryon vector form factor f 2 {q^) is in 
general suppressed by two more powers of the dibaryon invariant mass (f‘ than fi{q‘^) and since 
~ is large in S —> AcAcK decays, we expect the contributions from Xcc coupled to Ac Ac 

through be small and hence can be neglected in our calculation. The term in the 

axial-vector decay amplitude can also be dropped since qy{—g^^ + q^q'^/m\) ~ 0 due to the fact 
that m\ ~ q^ within the phase space. The decay amplitude Ax^^ then becomes 


_ _ Gf 

^v{A){B AjAc K ) = —^VcsVcba2mv(A)fv(A)Ui 




with 


Mv 

Ms 

Ma 


M p 


( hXXv \ 

2i^f^(g^) 2 2 ■ -’ 

\q^ - my + i myry j 

0 , 

2Ff^(g2) 


-h^AcKcA 


q^ - m\ AimA^A j ’ 


_ ,4.) ( 1 


x2mA. 


_hXAcA 


q^ - m\ XimA^A j 



(15) 


The minus signs in front of the strong couplings }iXXV,A from the minus sign of the part 
of the Xcc propagator in Eq. m- 
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FIG. 2: Branching fractions of —> A+A^, AT ) as a function of the mass and the decay 

width of the intermediate resonance Xcc when it is (a) a vector {Xcs = V) and (b) an axial- 


vector {Xcc = A) particle with /y • = 4 GeV. 


' AcAc-A 


We take /y • = fA ■ = 4 GeV and show in Fig. El the plots of branching ratios 

for each kind of resonance as functions of both the mass and the width of the resonance. The 
branching fractions depend on {fxcc ' and {fxcc ' proportionally. We notice that 

the width of the axial-vector resonance F^ is more than an order of magnitude smaller than Fy 
when both By and Ba are around the experimental value B ^ 7 x 10“^. This is due mainly 
to the smallness of Ba which suffers from the destructive interference between comparable A4 a 
and A4p contributions in the decay rate. Fig. 0 shows the decay rates from Ma, Xip and from 
ReAdy^Adp alone without taking into account the resonance effect. One can see that the Ma and 
Mp contributions are comparable while the interference term ReAd^Adp gives almost twice the 
negative contribution of either Ma or Mp, resulting in a large cancellation with Ma and Mp. 
Therefore, in order to counteract this cancellation, one needs a smaller width of the resonance such 
that |l/((7^ — m\-\-imATA)\~‘^ becomes more singular in the allowed range of q^. There is, however, 
no such interference found in By as My stands alone in the vector-induced decay amplitude after 
taking = 0 in Eq. (|T^ . 

Therefore, the above analysis seems to imply the existence of a narrow hidden charm bound state 
with a mass of order 4.6 ~ 4.7 GeV that couples strongly with the charmed baryon pair. Recall 
that many new charmonium-like resonances with masses around 4 GeV starting with X(3872) 
[l^ and so far ending with Y (4260) have been recently observed by BaBar and Belle. These 
charmonium-like states are above the DD threshold but below the two-charmed-baryon threshold. 
The new state we have put forward is just marginally above the AcAc threshold. In principle, this 
new state can be searched for in B decays and in pp collisions by studying the mass spectrum of 
or A,A,. 
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FIG. 3: Decay rates contributed from A4v (short-dash), Ma (solid), A4p (long-dash), sum of the 
previous two (dotted) and from 2ReM.A-M*p (dot-dashed) alone without taking into account the 
resonance effect, i.e. {q^ — m\ + iniA^ a)~^ is replaced by any constant in each term. 

III. TWO-BODY DECAYS 


The two-body doubly charmed baryonic B decays B~ —> and B^ —> receive contri¬ 

butions from the internal W-emission (see Fig. 0]) and weak annihilation. The latter contribution 
can be safely neglected as it is not only quark-mixing but also helicity suppressed. It should be 
stressed that, in contrast to the internal VF-emission in mesonic B decays, internal IF-emission in 
baryonic B decay is not necessarily color suppressed. This is because the baryon wave function is 
totally antisymmetric in color indices. One can see from Fig. |l]that there is no color suppression 
for the meson production. In the effective Hamiltonian approach, the relevant weak Hamiltonian 
is 


Weff = ^F,feC(ciOi + C2O2) ^ ^V,bV:,{ci - C2)Oi, 


(16) 


where Oi = {cb){sc) and O 2 = {cc){sb). In the above equation, we have used the fact that 
the operator Oi — O 2 is antisymmetric in color indices (more precisely, it is a color anti-triplet). 
Therefore, the Wilson coefficient for the tree-dominated internal kF-emission is ci — C 2 rather than 
a 2 = C 2 + ci/3. This is indeed the case found in the pole model calculation in |lf)l| . 

Since the internal VF-emission in Fig. I^is not factorizable, it is difficult to evaluate its amplitude 
directly. Pole model has been applied in [l^ to compute B — > BiB2- However, the strong coupling 
involved in this model is unknown and hence it has to be fixed from other processes, e.g. the 
3-body baryonic B decays. Since the CKM angles for B^ — > H^Aj and ^ same in 

magnitude (but opposite in sign), the pole model does not explicitly explain why the former has a 
rate much larger than the latter. In particular, the dynamical suppression of A^p relative to H^A“ 
is not clearly manifested in the pole model calculation. In order to understand why H^A“ 3> A^p, 
let us re-examine Fig. 0] 

There are several possibilities for the quark-antiquark pair creation in Fig. In one case, qq is 
picked up from the vacuum via the soft nonperturbative interactions so that it carries the vacuum 
quantum numbers ^Pq. It is also possible that the quark pair is created perturbatively via one 
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(a) 



c 

s 


d 


^0 



FIG. 4; 


and (b) i? —> as proceeding through internal VF-emission diagrams. 


gluon exchange with one-gluon quantum numbers ^Si. It is not clear which mechanism, the ^Pq or 
model, dominates the 2-body baryonic B decays, though in practice the ^Pq model is simpler. 
Since the energy release is relatively small in charmful baryonic B decay, the ^Pq model seems to 
be more relevant. In the present work, we also consider the possibility that the qq pair is produced 
via a light meson exchange. The qq pair created from soft nonperturbative interactions tends to 
be soft. For an energetic proton produced in 2-body B decays, the momentum fraction carried by 
its quark is large, ~ 0 ( 1 ), while for an energetic charmed baryon, its momentum is carried mostly 
by the charmed quark. As a consequence, the doubly charmed baryon state such as 'E^Ac has a 
configuration more favorable than Acp. 

In order to evaluate Fig. |3]for the decay B —> HcAc, we need to know the distribution amplitudes 
of the charmed baryon Be and the B meson. For the wave functions of Be = Ac, they have the 
forms 




^abc 


(Ac(p')|c“(4)4(^2)d^(4)|0) = 


he 


6 4 

/Ac r- 


^HO(P) C S5(?^ + mSc) ^E0izi,Z2,Z3), 
^ j a I J 7p ^ 


g 4 ['f^Ac(y)]„ [(/-"^Ac)75C']^^4’a^(7,4,4)> (17) 

where c, q and d are the quark fields, o, b, and c the color indices, a, (3 and 7 the spinor indices, C 
the charge conjugation matrix, and the decay constant. Following 1^ we can write 


'^Sc{zi,Z2,Z3) = y'[dx][d^fe_L] e*^*^''^'4'Hc(a;i,a:2,a:3,ki_L,k2A,k3A), 

4'Ac(z1,4>4) = [ [dx'][d'^k'2_]d^’^hiq,^^(^x[,X2,X3,k'i±,k'2±y3±), 


(18) 

where ,0±) is the momentum of HJ) (Ay) and ki\ k ^2 ^ ^ 3 ^ the momenta of the 

constituent quarks of the baryons, which are taken to be jlTi. l 22 || 


ki = (xip+,p ,ki_L), k[ = {p'^,x\p' ,k 3 _L), 

^ 2 ( 3 ) = (2^2(3)4''*', 0,1^2(3)^) , ^2(3) ~ (0) ^2(3)1^ >^2(3)±) 

[dx^'^\ = dxi^dx2^dx{^h ^1 — , 

[d‘^kx'\ ~ d‘^ k^^d"^ k2xd'^ k^]_6‘^ (^ij_ + 1 ^ 2 ! “I" l^sl) ’ 


(19) 


( 20 ) 















FIG. 5: Soft q{e)q{e) [= {dd), dd {uu) and du (ud)] as produced through the cr, vr^ and vr^ meson 

exchanges. Inside the parentheses are the momenta that the constituent quarks are carrying. Note 
the arrows on the quark lines do not represent the momentum flows. 


with Xi being the momentum fractions associated with the baryon, and kj_L the corresponding 
transverse momenta. Note that for simplicity, light quark masses are neglected in Eq. (Uni. The B 
meson wave function is expressed as 

£ pi 

(0|9t(a)6J(22)|B(p)> = -i^^|tf + mB)75W| < . (21) 

with P2 ^ ((1 -OPb. (1 - ^)Pb.O±). Pi ^Pi ^ (^PB.^Pfj.Op) andp^ 

The B~ —> decay amplitude now consists of three parts corresponding to the exchange of 

the a, TT^ and 7r“ between the soft qq quark pair and the spectator as shown in Fig. 

A{B ‘=‘cA(, ) = Afj + A^o + A^- , 


where, for instance, in the case that the decay proceeds through the exchange of a or 


C(,o, = - C2) 

x(=;K(0)4(0)J;;(t)|0)(A2|4(0)p:;,(/)<;;,(t)|0)(0|f4(t')i>?(0)|B-) 

^ [7/2(1 — 75)]q,p [7^(1 ~ 75)1/35 r^yFj^jj/ (22) 

d^ j[dx\[dx'][£‘k±\[(f‘k'^{2T:)^ 

a2 

^u{TT)qq 


Y* /s/Sc/Ac (Cl - C 2 ) 

^2 64 18 




+ r<xw 


XUS. 


T ^( 7 ( 71)75 




(23) 


with Fj^,^/ = r.yy = 1 for a and = —r.yy = i 75 for -k^ in Eq. (12211 . Note that the factor of 1/18 in 
Eq. (1231) is the color factor. For the exchange, the terms g^qq, r.yy and (Ac|c^(0)ny (z')(iy (2;)|0) 
in A^o are replaced by V^g^rqq, -rT,y and (Ac|c5(0)fiy (/)uy(2;)|0) = -(Ac|c5(0)uy(2;)iiy(2:')|0), 
respectively, where Eq. (El) has been used. Due to the symmetry property of the Ac wave function 
given in Eq. (P|l. the vr^ contribution is the same as the one except for an enhancement 
factor of (\/2)^ arising from isospin. Due to the tiny mass difference between vr^ and 7r“ and the 
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extremely narrow widths of these two particles, we have A^- j = 2 to a very good precision. 
The momentum labels of the quarks in He and Ac are depicted in Fig. 1^1 Pb and pi are the momenta 
of the b quark and the light spectator quark of the B meson defined after Eq m, respectively, 
9a{-K)qq is the Coupling of the a (vr) meson with the qq pair, and 


a„ = -4m3„ 


(ms + - ml^ 


b„ = -4m3„ 


(m-A, 




m| 


(24) 


for fj and 


a-jr = dms.. 


for both TT^ and tt . 


ryiB - (m-H, - 


6^ = dms.. 


ml - {ruB - mA,y 


(25) 


Our results are consistent with heavy quark effective theory. It has been shown that in the 
heavy quark limit, the decay amplitude can be expressed as u [A + v with [l^ 


A = 2^/ms [Q!(ri - r2) -/?], S =-2^/mA [a(ri + r2) -/?], (26) 


where ri = m-.^/mB, T 2 = mA^/ms and a, f3 are two unknown parameters. It is easily seen that, 
Eqs. (dU), (1^^ and agree with each other after setting a = + a-^r, (3 = (3^ + 13-,^, with 

Ufj = (3(j (X 1 + ri + r 2 for a and ajr = /Ijr oc 1 + ri — r 2 for both vr^ and 7r“, where the overall 
coefficients can be easily determined from Eqs. (ESI), (El and (El- 

To proceed with the numerical calculations, we need first deal with the delta functions that 
impose constraints on the integration limits as well as relations between integral variables. We 
show in the Appendix the decay amplitude and the integral variables as a result of the delta 
function integrations. The wave functions are adopted to be 


^BSxi,X 2 ,xy) = j [(fkA]^Byxi,kiA) = iVB,xiX 2 X 3 exp 


mt 


mn 


m% 


2(3^xi 2f3^X2 213^X2, 


with 




for the B meson 


2l| . The wave functions obey the normalization 


T]T(xi,X 2 ,X 3 ) = 1, J (ix4>s(x) = 1. 


(27) 


3^Bc 

k?L + ml 

ki L + ml 

+ ml 

(28) 

(27r/32)2 

and 

2/32x1 

2/32x2 

2/32x3 

$b(x) = Nbx‘^{1 — x)^ exp 

1 x^m^ 

2 uJb 


(29) 


(30) 


The other input parameters are specified as follows. The decay constants for the charmed 
baryons can be related to that of the A;, baryon via the relation 0 




(31) 


valid in the heavy quark limit. Using /a^ = 2.71 x 10“^ GeV^ obtained from a fit of the PQCD 
calculation for A;, —> Ac decays to —> AcW) 221, it is found that /a^ = 6.7 x 10“^ GeV^ and 
= 6.2 X 10“^ GeV^, which are roughly ~ 1.3 — 2.3 times that of the results from QGD sum 
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TABLE I: Predictions on the branching ratios of B~ —> and —> E+A“ decays. The first 

and second errors come from the theoretical uncertainties in the parameters (5 and ujb, respectively, 
which are taken to be /3 = 1.20 ± 0.05 GeV and cob = 0.40 ± 0.05 GeV, and the third error from the 
baryon decay constants. Results shown in second and third raws are from tt or a exchange alone, 
respectively. 


Mode 

Theory (10 

Expt (10-3) 

Mode 


Theory (10 3) 

Expt (10-3) 

1 0 

|<1 

0 0 

[I] 

T 

1 

\ r \ q-|-0.6-|-5.1-1-6.1 

) ^*"^-0.6-1.9-1.9 

« 4.8 

B{'^ - 

^ S^A-) 

pi-h0.5-|-4.7-|-5.6 
^•^-0.6-1.7-1.7 

Ri 1.2 

10 

i<i 

0 u 

[I] 

T 

1 

\ 10+0.5-1-4.2+5.0 

Jtt J--o_0.5_i.6-1,6 


B{B° - 


17+0.5-1-3.9+4.7 
'-0.5-1.4-1.5 


10 

i<i 

0 0 

[I] 

T 

1 

N n 9+0.0-I-0.4+0.6 

)(T U.Z_Q ^_Q 2 


B{^ - 

Jo- 

n 9+0.0-1-0.4+0.6 
^•^-0.0-0.1-0.2 


rules [^. ^ In 

our calculations we 

shall employ the value of the A;, baryon decay constant which 


is in the middle of the range that has the PQCD value as the lower bound and the highest value 
from QCD sum rules as the upper bound. The deviations of the lower and the upper bounds from 
this central value are then taken as one of the theoretical errors in our model. 

For the B meson, we use fs = 0.2 GeV. For the coupling 5o-(7r)gqj the linear sigma model leads 
to QaNN = V^’iTT'N/f-K with /jr = 132 MeV, and g^NN = V^'fnNgA/f-K with qa — 1.25 from the 
Goldberger-Treiman relation. Hence, it is reasonable to take ga{T^q = 9a{7T)NN/^ ™ the constituent 
quark model. For the a meson, we use To- ~ rria = 600 MeV |^. The constituent quark masses 
appearing in the Be wave function are taken to be rhu = = 0.33 GeV and mg = 0.55 GeV j^, 

while ific = is employed. 

The decay rate is given by 


KK) = P^f|4|2 {mB+mi+m2fpl , „|2 [(-^i + "^i)(-g2 + m 2 ) + | 

^ ^ 47rj {El + mi){E2 + m2)m% {Ei + mi){E2 + m2)m% J’ 

where Pc is the c.m. momentum, Ei and ruj are the energy and mass of the baryon Bi, respectively. 
The results of calculations are summarized in Tabled The theoretical errors come from the uncer¬ 
tainties in the parameters /3, uJb, which are taken to be /3 = 1.20±0.05 GeV LVb = 0.40±0.05 GeV, 
and the baryon decay constants. It is clear that the pion exchange gives the dominant contribution 
owing to its narrow width. The prediction is in agreement with experiment for H^A“, but a bit 
small for H[?A“ 

The above calculation is not applicable to the two-body decay ^ —>■ A+p with one charmed 
baryon in the final state. This is because two hard gluons are needed to produce an energetic 
antiproton: one hard gluon for kicking the spectator quark of the B meson to make it energetic and 




^ The Ab decay constant is found to be in the range (2.0 — 3.5) x 10“^ GeV^ in QCD sum rules I 23 . It 
differs from the Af, decay constant in this work by a factor of Af, mass. After normalizing it to having the 
same dimension as the decay constants in this work, the above range turns out to be about 1.3 — 2.3 times 
that of the depav constant from PQCD. 

^ As shown in |20j, the parameter (3 is of order 1 GeV for light baryons. Just as the meson case [llj, /3 
should become larger for the heavy baryons. 
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the other for producing the hard qq pair. The pQCD calculation for this decay will be much more 
involved (see e.g. for pQCD calculations of Af, —> AJ/if;) and is beyond the scope of the present 
work. Nevertheless, it is expected that r(i? —> BcN) <C r(i? —> HcAc) as the former is suppressed 
by order of af. This dynamical suppression effect for the A(^ production relative to He Ac has been 
neglected in the previous studies based on QCD sum rules 


]y| and on the diquark model |in| . 


IV. CONCLUSIONS 

In this work we have studied the two-body and three-body doubly charmed baryonic B decays, 
namely, B —> HcAc and B —> AcAcK, aiming to resolve the puzzles associated with them. We 
point out that the suppression from the A'^A~K~ phase space is so strong that A+A“ pair is 
unlikely to be produced dominantly through the direct three-body decay processes. Nevertheless, 
the great suppression for the A'^A~K production can be alleviated provided that there exists a 
narrow hidden charm bound state with a mass near the AcAc threshold, of order 4.6 ~ 4.7 GeV. 
This new state that couples strongly to the charmed baryon pair can be searched for in B decays 
and in pp collisions by studying the mass spectrum of or Ac Ac. 

The doubly charmful decay such as B ^ HcAc has a conhguration more favorable than the singly 
charmful one such as B^ —> Aep even though they have the same CKM angles in magnitude. This is 
because no hard gluon is needed to produce the energetic HcAc pair in the former decay, while two 
hard gluons are needed for the latter process. Therefore, Aep is suppressed relative to HcAc due to 
a dynamical suppression from C>(af). Assuming that a soft qq quark pair is produced through the 
a and tt meson exchanges in the conhguration for B —> HcAc, it is found that B{B HcAc) ~ 10“^. 
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Note added: After this paper was submitted for publication, Belle published the updated version 
of in which the spectrum of the AcAc pair in the B AcAcK decays is shown for the hrst time. 
According to Belle’s observation no new resonance with a mass near the AcAc threshold was found 
(see Fig. 3 in version 2 of Q]). This implies the failure of naive factorization for this decay mode 
and may hint at the importance of nonfactorizable contributions such as hnal-state effects. For 
example, the weak decay B —> followed by the rescattering —> A^A^K Q or the 

decay B —> HcAc followed by HcAc —> AcAcK may explain the large rate observed for B —> AcAcK. 
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APPENDIX A: 


After integrating out the delta functions in Eq. (US}, we obtain, taking decay via a as an example, 

A^{B- ^ hOA-) 


= - C2)(27r)42 


rp' /pb /■i-«Ps/p+ dx 2 [^Pb/p' dxi 


di 


Jo 


p 


2 

J- 


x( 27 r) 


dk 


■ 3 ± 




r2-K 




d^23 

Jo 

(do-gq) 


dk 


.12 

■ 3 ± 


r2-K 


de 


33 ' 


■ U'ESaa + ba'y5)vAc , 


(Al) 


pI- ml + i 

where ,0±) is the momentum of (A“) and ki \ ^ k^ the momenta of the 

constituent quarks of the baryons, which are taken to be 

ki = (xip+,p",ki_L), k[ = (y+,xiy",ki_L), 

^2(3) = (2^2(3)P^) 0) k2(3)_i_) , ^2(3) = (0>^2(3)^'^ >^^2(3)±)- 

The integrations of the delta functions yield the following relations 

ipt 


J) 1 V) 

X' ' I ‘ 

1 - X U/ O db 


(>) A') 


3 ’ 


3^3 = 


/ iPB 
H- 


P 


J- 


^2 ) 


kS'i = - (k?i + kW ) , k'A = - (ksx + k'^ k;^ = ksA 

and the limits of integrations as shown in Eq. (EB). Thus, 

^^Pb 


(A3) 


Pa — (^3 + ^ 3 ) — 


p-i 


p 


— x'A p'^p — ^k|_|_ + k 3 \ + 2 |k 3 _L| |k 3 _|_| cos ^ 33 '^ , 


k?± = k2_L+k|_L+ 2|k2A||k3A|cos6'23, 
k2± = k|_L + kfj_ + 2|k3_L||k(5_L|cos(933/, 


(A4) 


where 023 and ^ 33 / are the angles of k 2 _L and of kg^ as measured against k 3 _i_, respectively. Note 
that the integration ranges of X 2 and X 2 are constrained by the delta function in Eq. (El, while 
the integration range of ^ is restricted by (1 — Cp~b/p^) — ^P~b/p'~ < 1 - 
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